UNCLASSIFIED 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
ELECTRONICS  RESEARCH  LABORATORY 


AR-OOI -667 


TECHNICAL  REPORT 
ERL-0077-TR 


THE  MEASUREMENT 
ON  THE  LIQUID 


OF  TEMPERATURE  INSTABILITIES  AND  THEIR  EFFECT 
PHASE  EPITAXIAL  GROWTH  OF  GALLIUM  ARSENIDE 


M.  Folkard 


SUMMARY 

Epitaxial  layers  grown  on  gallium  arsenide  substrates  using 
the  steady  state  temperature  gradient  form  of  liquid  phase 
epitaxy  have  a  much  larger  growth  rate  than  that  calculated 
for  diffusion-limited  conditions.  In  addition,  and  of  major 
importance  for  device  fabrication,  there  is  a  large  variation 
in  epilayer  thickness  over  the  surface  of  a  single  sample,  and 
non-reproducibility  from  sample  to  sample. 

The  present  report  describes  various  measurements  of  the 
temperature  and  its  local  variation  within  the  gallium  melt 
solution,  and  attempts  to  relate  the  observed  variations  and 
instabilities  to  the  epitaxial  layers  actually  grown. 
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1  .  INTRODUCTION 


Cooling  a  gallium  melt  which  is  saturated  with  arsenic  will  result  in  the 
epitaxial  deposition  of  a  monocrystalline  layer  onto  a  single  crystal  gallium 
arsenide  substrate  placed  in  contact  with  the  melt.  This  occurs,  in  part, 
because  the  solubility  of  a  dilute  constituent  in  a  liquid  solvent  decreases 
with  decreasing  temperature.  A  variation  of  this  method,  which  allows  epitaxial 
growth  to  take  place  at  a  constant  temperature,  uses  a  temperature  differential 
imposed  across  the  melt  solution  by  a  split  furnace  whose  upper  element  is 
maintained  at  a  lower  temperature  than  the  bottom  element.  This  temperature 
gradient  across  the  melt  provides  a  driving  force  for  the  diffusion  of  arsenic, 
dissolved  from  a  gallium  arsenide  source,  through  the  melt  to  the  cooler  sub¬ 
strate  where  it  is  deposited  as  an  epitaxial  layer  of  gallium  arsenide. 

The  growth  rates  predicted  under  specified  conditions  using  a  diffusion 
model  are  much  smaller  than  the  rate  of  epilayer  growth  which  is  experimentally 
observed(ref . 1 ) .  In  addition,  and  of  major  importance  for  device  fabrication, 
there  is  a  large  variation  in  epilayer  thickness  over  the  surface  of  a  single 
sample,  and  non-reproducability  from  sample  to  sample. 

The  present  report  describes  various  measurements  of  the  temperature  and  its 
local  variation  within  the  melt  solution,  and  attempts  to  relate  the  observed 
variation  in  epitaxial  layer  heights  to  these  temperature  variations. 


2.  THE  CALCULATION  OP  GROWTH  RATES 

A  mathematical  formulation  for  the  growth  of  epitaxial  layers  by  the  steady 
state  method  may  be  derived  by  assuming  that  diffusion  of  gallium  arsenide 
through  the  melt  is  the  rate-limiting  process. 

Such  a  treatment(ref .1 )  predicts  a  growth  rate: 
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temperature  gradient  within  the  melt  solution 

temperature  of  the  source  material 
temperature  of  the  substrate  surface. 
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To  estimate  the  rate  of  epilayer  growth  it  is  necessary  to  obtain  relevant 
values  for  the  slope  of  the  liquidus  curve,  the  diffusion  coefficient,  and  the 
solution  temperature  gradient. 

2.1  Slope  of  the  liquidus  curve 

dX^  ( T ) 

The  slope  As _  7  of  the  liquidus  curve  is  obtained  from  the  solubility 
dT 

data  of  Hall(ref.2)  which  receives  independant  verification  from  Hsieh(ref.3) 
in  the  limited  temperature  region  near  800  C.  A  smooth  curve  drawn  through 
Hall's  experimental  data,  with  almost  every  experimental  data  point  lying 
actually  on  the  smooth  curve,  yields  the  values 

T  =  700°C,  liquidus  slope  =  0.90  x  10  ^/°C 
T  =  800°C,  liquidus  slope  -  2.78  x  10  ^/°C 


(These  values  are  considered  accurate  to  within  a  few  percent.) 

2.2  Diffusion  coefficient  D 

There  have  been  no  direct  measurements  of  the  diffusion  of  As  or  GaAs  in 
gallium  solution.  Estimates  for  D  have  been  made  by  comparing  experimental 
growth  rates  with  the  predictions  of  various  mathematical  theories  based  on 
diffusion-limited  growth,  and  obtaining  "best  fit"  values  for  D.  There  is 

_5 

a  large  spread  in  the  values  so  obtained,  for  example  from  1x10  to 

10  x  10“5  cm2/s  at  800°C,  with  some  investigators  claiming  no  temperature 
dependence  and  others  claiming  an  exponential  dependence  (ref. 3  to  7). 

The  values 

T  =  700°C  D  =  3  x  10-5  cm2/s 

T  =  800°C  D  =  5  x  10~5  cm2/s 

are  consistent  with  the  general  trend  of  the  majority  of  these  references, 
and  so  have  been  used  in  the  present  case. 

AT 

2.3  Solution  temperature  gradient  _ s 

Ax 

To  determine  the  temperature  gradient  within  the  melt,  a  pyrolytic  graphite 
probe  was  constructed  containing  a  chrome 1-alumel  thermocouple.  Full  details 
of  the  probe  and  the  methods  of  measurement  are  given  in  Section  4.  The 
repeatability  of  these  measurements  was  poor,  as  can  be  seen  in  Tables  1  to 
3,  but  they  did  indicate  that  the  solution  temperature  gradient  was  relatively 
independent  of  the  actual  furnace  temperatures  for  a  given  differential  AT 
between  top  and  bottom  furnace  elements.  For  an  external  differential  of 
50°G,  the  solution  temperature  gradient  was  found  to  be  approximately  4.2 
C/cm  at  700°C  (Table  4J  and  4.6  C/cm  at  800°C  (Tables  1  to  3). 

2.4  Predicted  growth  rates 

Using  the  values  shown  below,  the  predicted  growth  rates  become 


1 .0  p/hour  at  700°C,  and 
5.3  p/hour  at  800°C. 
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.  t(°c) 

D  (cm^/s) 

(slope  x\s  (T)  .  t) 

sub 

AT 

(  Ax  )  C/cm 

700 

3  x  10-5 

0.90  x  10“4 

4.2 

800 

5  x  10“5 

2.78  x  10"4 

4.6 

3.  OBSERVED  GROWTH  RATES 


The  following  results  apply  for  a  steady  state  temperature  differential  AT 
of  50°C,  and  indicate  the  order  of  magnitude  of  the  growth  rates  obtained  for 
2  hours  at  800UC  or  4  hours  at  700  G.  Before  presenting  the  observed  growth 
rates,  some  details  on  the  apparatus  and  measuring  technique  are  relevant. 

3.1  The  growth  vessel 

The  growth  vessel  of  vitreous  carbon  is  shown  in  figure  1 ,  where  it  can 
be  seen  that  the  source  block  is  held  down  by  a  graphite  peg  of  length  31  mm 
and  square  cross  section  of  side  6  mm.  The  substrate  (5  mm  x  10  mm)  is  held 
by  suction  onto  a  pyrolytic  graphite  holder  of  length  36  mm  and  area  13  mm  x 
7  mm  at  substrate  and  1 9  mm  x  9  mm  at  top.  Movement  of  the  substrate  and 
holder  is  manually  controlled  from  outside  the  furnace  by  a  long  quartz  rod, 
and  growth  is  initiated  (and  terminated)  by  rotating  the  substrate  holder 
down  into  (or  up  out  of)  the  melt  solution. 

3.2  Measurement  of  the  epilayer  heights 

Spot  heights  of  the  grown  layers  were  found  by  breaking  each  sample  into 

sections  along  ( 1 1 0 )  planes  (shown  in  figure  2  by  faint  lines)  then  staining 

for  1 5  to  40  s  with  a  solution  of  HF  :  HN0„  :  H_0  =  1  :3:4.  Differences  in 

o  2 

electron  affinities  in  the  differently  doped  regions  affect  the  rate  at  which 
GaAs  reacts  with  the  staining  solution,  and  the  difference  in  etching  rates 
causes  the  junction  to  become  visible.  A  Zeiss  KK08  reflected  light  micro¬ 
scope  was  used  with  a  Nomarski  differential  interference  contrast  attachment 
and  a  magnification  of  1 28X.  Calibration  of  the  eyepiece  scale  yielded  a 
value  of  6.8  y  per  smallest  division  at  this  magnification,  so  that  the 
resolution  of  the  tabulated  layer  thickness  is  within  a  few  microns  in  each 
case. 

3.3  Observed  epilayer  heights 

Figure  2  is  a  schematic  representation  of  the  epilayer  thickness  for 
eight  individual  samples.  The  left  hand  side  of  each  diagram  represents 
the  edge  of  the  substrate  which  was  closest  to  the  graphite  peg  used  to 
hold  the  solute  block  in  position  in  the  bottom  of  the  crucible.  Although 
the  results  from  each  sample  were  different  the  general  trend  shows  a 
thicker  grown  layer  on  the  side  of  the  substrate  nearer  the  graphite  peg. 

This  suggests  that  the  portion  of  the  melt  near  the  graphite  peg  may  be 
cooler  than  elsewhere.  The  variable  epilayer  heights  over  a  single  sample, 
and  the  non-repeatability  from  sample  to  sample  quite  strongly  suggest  that 
processes  other  than  simple  diffusion  of  gallium  arsenide  through  the  melt 
are  taking  place.  It  is  also  evident  from  figure  2  that  the  observed  growth 
rates  are  greatly  in  excess  of  the  values  calculated  in  Section  2.  A  very 
rough  estimate  for  the  "average"  growth  rate  for  the  samples  shown  in  figure 
2  would  be  5  y/hour  and  15  y/hour,  at  700  C  and  800  C,  compared  to  the 
calculated  values  of  1  y/hour  and  5  y/hour  respectively.  However,  the 
dominant  feature  is  the  non-uniformity  of  the  growth. 
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3.4  Discussion 

Uniform  growth  rates  much  greater  than  those  predicted  from  the  diffusion 
model,  and  more  in  accord  with  the  "average"  magnitudes  observed,  might  be 
explained  by  considering  either 

AT 

(a)  a  temperature  gradient  _ _  of  about  15  C/cm  within  the  solution, 

Ax 

or  (b)  a  much  greater  value  for  D  than  that  used  in  the  calculations  of 
Section  2. 

Possibility  (a)  is  improbable  because  the  total  temperature  differential 
of  50°C  occurs  over  a  distance  in  excess  of  6  cm,  and  measurements  later 
described  in  the  present  report  (Table  5)  indicate  the  temperature  gradient 
in  the  hydrogen  atmosphere  to  be  8.0°C/cm.  Possibility  (b)  is  also  unable 
to  account  for  the  increased  growth  rate.  Even  considering  the  largest 
value  of  10  x  10“  ^  cm^/s  estimated  for  D  from  the  wide  range  of  values  in 
references  3  to  7,  the  predicted  growth  rate  using  the  diffusion  model  is 
still  only  ~1 0. 6  y/hour. 

However,  the  dominant  feature  observed  from  the  various  samples  in 
figure  2  is  the  non-uniformity  of  grown  layer  heights  over  the  area  of  each 
sample,  which  can  in  no  way  be  accounted  for  by  the  simple  diffusion  model. 

A  plausible  explanation  for  the  observed  situation  is  that  kinetic  processes 
such  as  convection  might  be  occuring  within  the  melt  in  addition  to  the 
diffusion  of  arsenic  through  the  melt.  To  test  this  hypothesis,  some 
variation  to  the  standard  experimental  procedures  of  Runs  1  to  1 5  were  made. 
In  these  earlier  experiments,  the  epitaxial  layers  had  been  grown  from  melts 
containing  30  g  of  gallium,  which  filled  the  crucible  to  7  mm  above  the 
solute  block,  while  the  maximum  height  of  the  substrate  above  the  solute 
block  was  5  mm.  In  Runs  16  to  20,  the  melt  was  increased  to  60  g  of  gallium, 
which  brings  the  liquid  surface  about  13  mm  above  the  solute  block,  so  that 
the  liquid  surface  is  now  at  least  8  mm  above  the  growth  plane.  Consequently 
there  is  a  large  volume  of  melt  above  the  substrate  and  we  might  expect 
convective  effects  to  be  increased  by  this  larger  melt  volume.  In  these 
circumstances  the  effect  on  a  substrate  at  700°C  for  4  hours  was  to  grow  a 
layer  on  the  side  of  the  substrate  near  the  graphite  peg  and  to  etch  material 
away  from  the  side  further  into  the  melt  (see  figure  3).  A  similar  result 
was  obtained  at  800°C  for  2  hours,  while  at  the  higher  temperature  a  "growth" 
time  of  4  hours  resulted  in  the  substrate  being  etched  all  over,  though  to  a 
lesser  extent  near  the  peg  than  remote  from  it.  The  substrate  -  nominally 
at  a  cooler  temperature  than  the  solute  block  -  had  been  partially  dissolved 
although  there  was  a  large  excess  of  solute  present.  In  addition,  the 
solute  block  exhibited  a  definite  pattern  of  channels  on  its  surface.  Par 
from  the  graphite  peg,  the  solute  block  had  been  completely  dissolved  while 
no  dissolution  of  material  whatever  had  taken  place  actually  under  the  peg. 
This  provides  definite  evidence  for  non-uniformities  within  the  growth 
environment,  and  makes  it  imperative  to  study  the  growth  kinetics  more 
closely.  The  first  step  in  determining  how  convection  and  associated 
processes  may  be  developing  is  to  determine  the  variation  in  temperature 
and  temperature  gradient  throughout  the  actual  melt.  This  involves  precise 
temperature  profiling  under  specified  conditions. 
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4.  TEMPERATURE  PROFILE  IN  THE  MELT 

Preliminary  measurements  using  a  single  thermocouple  probe  are  first  described, 
followed  by  the  more  detailed  measurements  in  which  a  differential  thermocouple 
probe  was  used.  In  each  case  the  monitoring  system  simulated  the  experimental 
situation  as  closely  as  possible,  aiming  to  determine  the  temperature  conditions 
during  an  actual  growth  run  rather  than  the  conditions  in  an  'unperturbed  melt. 

Since  the  thermal  conductivity  of  the  pyrolytic  graphite  substrate  holder  itself 
must  alter  the  temperature  gradient  within  the  solution  (the  electrical  and 
thermal  properties  of  pyrolytic  graphite  are  markedly  anistropic)  the  probe  was 
made  of  the  same  dimensions  as  the  substrate  holder,  and  the  same  quartz  dipping 
rod  arrangement  was  used. 

4.1  Single  Thermocouple  Probe 

The  initial  temperature  measurements  were  taken  with  a  chrome 1-alumel 
thermocouple  of  0.025  inches  diameter  wire.  Quartz  capillary  tubes  were 
placed  over  each  of  the  leads  to  prevent  contact  with  each  other  and  the 
graphite  probe.  The  "hot"  junction  of  the  thermocouple  was  encased  in 
quartz  which  was  ground  flat  and  held  against  the  bottom  of  the  graphite 
probe  by  a  small  pyrolytic  graphite  screw,  in  the  region  where  the  substrate 
would  normally  be  positioned.  Although  this  screw  prevented  the  thermo¬ 
couple  from  reaching  closer  than  2  mm  to  the  crucible  base,  it  was  felt  a 
not-too-unreasonable  perturbation  for  a  first  attempt,  since  actual  growths 
are  carried  out  at  distances  greater  than  2  mm  above  the  crucible  bottom. 

The  probe  was  dipped  into  a  Ga  melt  of  depth  about  18  mm  and  raised  or 
lowered  to  simulate  the  geometry  used  for  epitaxial  growth.  The  melt 
temperature  T  (x)  was  monitored  as  a  function  of  distance  x  above  the 

bottom  of  the  crucible  using  a  Rikadenki  B-24  recorder  which  was  calibrated 

with  a  Cambridge  Model  44228  potentiometer  bridge.  Comparison  of  T  (x)  at 

s 

different  positions  within  the  melt  allowed  a  temperature  profile  to  be 
obtained. 

To  increase  the  accuracy  of  measurement,  an  off-set  voltage  arrangement 
using  the  Cambridge  potentiometer  was  set  up  as  shown  in  figure  4.  It  was 
possible  by  this  means  to  make  differential  temperature  measurement  with 
great  accuracy  because  only  short  term  stability  of  the  thermocouple  and 
that  of  the  off-set  voltage  are  of  first  order  significance  in  determining 
the  accuracy  of  the  difference  between  any  two  of  the  readings. 

Although  each  individual  set  of  measurements  was  reasonably  self-consistant, 
the  repeatability  of  measurements  from  day  to  day,  or  even  in  successive 
raisings  and  lowerings  of  the  probe,  was  poor.  The  readings  from  a  typical 
experiment  are  shown  in  figure  5,  while  Table  1  demonstrates  the  variation 
in  results  observed  from  day  to  day.  For  the  conditions  under  which  most 
of  the  epitaxial  growths  were  carried  out  (T^  =  769°C,  bottom  =  818°C, 

AT  nominally  50°C)  the  temperature  gradient  using  this  probe  was  found  to 
be  of  the  order  of  1  °C/ cm. 

In  these  experiments,  no  correction  was  applied  to  the  results  to  allow 
for  the  change  in  relative  depth  of  solution  with  displacement  by  the 
measuring  probe  (in  any  case  corrections  would  be  of  minor  significance  here 
where  a  large  volume  of  Ga  was  used,  although  it  might  be  more  significant 
in  the  situation  where  only  a  small  volume  of  melt  was  used).  Of  more 
importance  however,  is  the  assumption  that  the  temperature  profile  is  not 
modified  by  moving  the  probe. 


ERL-0077-TR 


-  6  - 


4.2  Differential  Thermocouple  Probe 


To  overcome  some  of  the  limitations  of  the  single  thermocouple  measure¬ 
ments,  a  differential  thermocouple  was  constructed  and  used  with  the  same 
graphite  probe  and  dipping  rod  arrangement  as  has  already  been  described. 

This  time  the  wires  were  of  0.005"  diameter,  which  caused  some  difficulties 
as  they  became  brittle  after  sustained  operation  at  high  temperature  in  the 
reducing  atmosphere  and  sometimes  broke.  When  this  occurred,  the  furnance 
had  to  be  cooled  and  the  entire  apparatus  dismantled,  the  wires  repaired  and 
the  apparatus  reassembled,  causing  considerable  delay.  The  vertical 
separation  of  the  two  "hot"  junctions  was  2.0  mm,  formed  by  a  length  of 
chromel  wire,  with  the  leadout  wires  made  from  alumel.  Again  the  thermo¬ 
couple  was  encased  in  thin  quartz  to  prevent  melting  by  the  hot  gallium 
solution,  and  the  leadout  wires  were  threaded  with  quartz  capillary  to 
prevent  contact  and  possible  shorting  with  each  other  and  with  the  graphite 
probe.  A  schematic  diagram  of  the  probe  assembly  is  shown  in  figure  6. 

To  increase  the  accuracy  of  measurements,  the  arrangement  shown  in  figure 
7  was  used.  The  Hewlett-Packard  41 9A  DC  null  voltmeter  output  terminal  was 
adjusted  to  read  (l.OOOO  i  0.0005)  volts  for  full-scale  deflection  on  any 
input  range  of  the  digital  voltmeter.  At  the  same  time  the  zero  and  gain  of 
the  Rikadenki  recorder  were  adjusted  on  the  1  volt  input  scale.  Once  set 
up,  these  instruments  were  found  to  remain  exceptionally  stable  over  several 
weeks  operation.  By  using  the  100  pV  range  on  the  HP  voltmeter,  the  25  cm 
scale  on  the  Rikadenki  recorder  then  corresponds  to  only  2.4°C,  so  that  this 
arrangement  allowed  the  detection  of  very  small  temperature  differences. 

Initially  the  probe  was  positioned  near,  but  not  actually  touching,  the 
bottom  of  the  crucible  in  a  melt  of  30  gm  Ga.  This  melt  fills  the  crucible 
(with  graphite  peg  present)  to  a  height  of  7  mm,  which  corresponds  to  the 
conditions  of  Runs  1  to  1 5  earlier  described.  With  the  probe  in  this 
position  the  temperature  differential  ATp  across  the  differential  thermo¬ 


couple  was  monitored  as  the  furnance  differential  AT  was  altered,  keeping 
the  "average"  furnance  temperature  at  793°C  (33  mV  on  the  chrome 1-alumel 

AT  AT 

monitoring  thermocouple).  Using  the  equation — -  0r 

Ax 

/AT 

relationship  (see  figure  8)  between  (  — -  )  and  AT  was  found  to  be 


0.2 


^/cm,  the 


Ax 


AT 


Ax 


=  3.85  x 


10-2 


AT 

It  should  be  observed  from  figure  8  that  although  linearity  is  shown, 

AT  \ 
s  \ 

-  j  and  AT  does  not  occur.  A  simple  treat- 

Ax  ' 

ment  by  McRae(ref.8)  which  considers  that  heat  transfer  from  the  furnace  to 
the  melt  takes  place  only  by  radiation,  predicts  that 


direct  proportionality  of 


AT  + 


terms  in  (  AT) 


2 


Rote  that  as  AT  goes  to  zero,  this  equation  predicts  that 


AT 
_ s 

Ax 


will 


also  vanish.  Figure  8  shows  that  experiment  does  not  support  this  prediction 
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and  we  are  forced  to  conclude  that  (  _s_  )  and  AT  are  not  directly  propor- 

Ax  ' 

tional. 

Moreover,  it  soon  became  apparent  from  these  measurements,  the  value  of 
AT  \ 

_ s_  )  when  AT  =  0  is  far  from  constant,  even  when  the  probe  is  not  moved 

Ax  / 

at  all.  In  fact  when  the  probe  differential  temperature  ATp  was  monitored 


over  extended  periods  with  AT  =  0,  a  linear  "drift"  with  time  was  found, 
of  the  order  of  several  microvolts/hour.  Furthermore  when  similar  measure¬ 
ments  are  carried  out  under  sensibly  constant  conditions  (probe  position 
unaltered)  over  periods  of  several  days,  the  observed  drift  may  be  positive 
or  negative.  Figure  9  shows  the  magnitude  of  the  baseline  drift  on  one 
occasion  over  a  1 2  hour  period  when  both  furnace  controls  were  set  to  33  mV 
(793°C).  From  Table  2,  where  measurements  were  made  over  a  5  day  period,  it 
can  be  seen  that  the  baseline  drift  varies  in  an  irregular  manner,  which 
makes  one  suspect  that  it  might  be  an  instrumentally  induced  effect. 

For  the  Series  A  measurements  recorded  in  Table  2,  the  temperature  con¬ 
troller  of  the  bottom  furnace  element  was  initially  set  at  34  mV  (818°C) 
while  the  top  controller  was  set  to  32  mV  (769°C),  producing  a  nominal 
temperature  difference  AT  of  50°C.  (These  conditions,  when  used  for 
epitaxial  layer  growths,  are  referred  to  in  the  present  report  as  "growth 
at  800°C".)  The  voltage  reading  (Tp)  ,  from  the  differential  probe  was 


recorded  under  these  conditions,  after  which  both  top  and  bottom  furnace 
controllers  were  reset  to  33  mV  (corresponding  to  793°C).  When  the  furnace 
had  stabilised  at  these  new  conditions  with  AT  =  0  another  probe  voltage 
reading  (Tp)^  was  recorded. 

The  voltage  from  the  probe  was  monitored  continuously  on  a  recorder,  as 
shown  for  a  typical  experiment  in  figure  10.  The  slope  of  the  trace  under 
constant  applied  conditions  allowed  an  estimation  of  the  baseline  drift 
during  the  time  between  the  measurement  of  (Tp)^  and  (Tp)^,  whilst  the 

difference  between  (Tp)^  and  (Tp)p  was  considered  as  the  real  measure  of 

ATp  induced  by  the  external  furnace  differential  of  AT  =  50°C.  (The 


voltage  overshoot  which  is  observed  from  the  probe  whenever  the  furnace 
controllers  are  altered  can  be  directly  explained  by  the  furnace  temperature 
overshoots  which  are  clearly  demonstrated  in  figure  11.)  The  value  of  ATp 


so  obtained  and  corrected  for  baseline  drift  was  then  converted  to  an  equiva¬ 
lent  temperature  gradient  in  °C/cm.  It  is  of  interest  to  note  that  the 
measured  temperature  gradient  increased  from  a  value  of  about  2.0  C/cm  during 
the  early  period  of  operation  at  800°C,  to  approximately  (4.6  1  0.4)°C/cm 
after  an  extended  period.  Additional  measurements  showed,  as  expected  from 
the  form  of  the  grown  layers  described  in  previous  sections,  that  the 
temperature  gradient  decreased  with  distance  away  from  the  graphite  peg. 

Another  set  of  measurements,  referred  to  as  Series  B  and  shown  in  Table  3 
were  carried  out  after  the  thermocouple  had  broken  and  been  repaired. 

Although  the  "relative"  temperature  gradient  for  the  first  day's  operation 
in  Series  B,  did  not  differ  appreciably  from  the  first  day's  operation  in 
Series  A,  and  the  magnitudes  of  the  baseline  drift  were  similar,  the 
"absolute"  probe  readings  were  very  much  greater  in  the  Series  B  measurements. 
The  probe  was  then  dismantled,  both  alumel  leadout  wires  were  sheathed  in 
quartz  capillary  and  threaded  down  a  cylindrical  copper  shielding  cable, 
which  was  in  turn  threaded  through  the  hollow  quartz  dipping  rod.  By  using 
this  careful  shielding  it  was  hoped  to 


(i)  reduce  the  noise  on  the  output  signal  of  the  differential 
thermocouple 
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(ii)  reduce  or  eliminate  the  variable  baseline  drift  on  the  output 

signal  which  occurred  under  supposedly  steady  furnace  conditions 

(iii)  reduce  or  even  eliminate  whatever  effects  were  causing  the  variable 
(and  sometimes  very  large)  readings  (T  on  the  differential 

thermocouple  probe  when  AT  =  0,  i.e.  when  no  temperature  differential 
is  applied  across  the  furnace. 

Although  the  shielding  reduced  the  noise  level  a  great  deal,  and  also 
eliminated  the  baseline  drift  once  the  probe  had  been  in  the  melt  for 
several  days  (Tables  4  to  6),  the  results  of  Series  C  measurements  did  not 
allow  any  quantitive  conclusions  to  be  made  regarding  the  effects  causing 
a  non-zero  probe  reading  (Tp)^  for  AT  =  0. 

All  the  measurements  in  Tables  2  to  4  were  made  at  a  nominal  "average" 
temperature  of  800°C,  with  either  the  top  furnace  controller  set  to  34  mV 
(818°C)  and  the  bottom  set  to  32  mV  (769°C),  or  the  top  and  bottom  set 
together  to  33  mV  (793°C).  Further  measurements  made  during  Series  C  with 
bottom  and  top  furnace  settings  of  30  mV  (721 °C)  and  28  mV  (673°G)  respectively 
or  both  elements  at  29  mV  (697°C)  are  presented  in  Table  5.  These  readings 
at  a  nominal  "average"  temperature  of  700°C  produced  a  value  for  the  solution 

temperature  gradient  ^s  of  approximately  4.2°C/cm.  The  results  of  4.6°C/cm 

Ax 

at  an  average  temperature  of  800°C  and  4.2°C/cm  at  an  average  temperature  of 
700°C,  show  that  (  ^~s  j varies  only  slowly  with  the  average  furnace  temperature 

V  Ax/ 

provided  that  AT  is  kept  constant. 

The  temperature  gradient  in  the  flowing  hydrogen  stream  above  the  crucible 
was  measured  by  raising  the  probe  vertically  out  of  the  melt  as  far  as  it  would 
go  towards  the  top  of  the  quartz  furnace  tube.  It  was  visibly  checked  that 
no  Ga  remained  on  the  probe  before  a  series  of  measurements  were  made  at  an 
average  furnace  temperature  of  800°G  and  a  differential  of  50°C,  in  similar 
fashion  to  those  made  when  the  probe  is  in  the  melt.  The  observed  value 
AT 

a2  of  8.0°C/ cm  shown  in  Table  6,  agrees  well  with  the  known  applied  AT 

Ax 

of  50°C  applied  over  a  6  cm  region,  and  so  provides  confidence  in  the  results 
derived  using  this  probe. 

However,  on  reviewing  all  the  information  shown  in  Tables  2  to  4,  it  is 
evident  that  in  each  case  the  early  measurements  yield  a  value  for  the 
solution  temperature  gradient  of  about  2°C/cm,  while  those  measurements 
taken  after  the  furnace  had  been  on  for  an  extended  period  (several  days) 
all  yield  values  close  to  4.6°C/cm.  Insufficient  time  remained  to  investi¬ 
gate  this  phenomenon,  but  it  possibly  is  consistent  with  the  observation 
that  about  24  hours  "equilibration"  time  between  the  melt  and  the  source 
material  is  necessary  before  good  epitaxial  layers  can  be  grown  routinely 
(ref.l). 

The  conclusions  drawn  from  these  temperature  probing  experiments  using 
the  differential  thermocouple  probe  are  as  follows: 


"Average"  Temperature 
T(°C) 

Furnace  temperature 
differential 

AT  (°C) 

Solution 

temperature 

Gradient 

AT  /Op/  \ 
s  (.  C/cm) 

Ax 

Temperature  gradient 
in  Hydrogen  atmosphere 

(  j)  C/0m 

700°C 

50°C 

4.2 

800°  C 

50°C 

4.6 

8.0°  C/cm 
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It  must  also  be  remembered  that (  _ s  ]  varies  within  the  actual  solution, 

'  Ax' 

tending  to  decrease  near  the  graphite  peg,  and  also  to  decrease  towards  the 
bottom  of  the  crucible  -  so  that  the  values  tabulated  above  are  not  of  high 
absolute  accuracy. 


4.3  Furnace  temperature  uniformity 

/  AT  \ 

Because  of  the  substantial  variation  in  the  value  for  I  _ s  J  obtained 

'  Ax/ 

from  day  to  day  with  the  single  thermocouple  probe,  some  measurements  of  the 
axial  variation  of  furnace  temperature  itself  were  carried  out.  Top  and 
bottom  sections  of  the  Marshal  split  furnace  cars  be  controlled  independently, 
and  each  section  consists  of  3  appropriately  connected  elements.  Temperature 
control  in  each  section  is  maintained  by  feedback  from  a  chrome 1-alumel 
thermocouple  fixed  near  the  mid-position  of  each  centre  element. 

The  temperature  variations  along  both  top  and  bottom  sections  of  the 
furnace  were  determined  by  moving  independent  chrome 1-alumel  monitoring 
thermocouples  axially  along  each  section,  and  measuring  the  output  voltage 
with  the  Cambridge  potentiometer.  As  shown  in  figure  12,  even  when  both 
controllers  are  set  to  the  same  temperature  of  793  C  (33  mV),  the  axial 
variation  is  significant,  with  a  steady  increase  in  temperature  out  to  a 
distance  of  about  12  cm  from  the  centre,  followed  by  a  rapid  decrease  as  the 
open  ends  of  the  furnace  exert  increasing  influence.  Although  the  bottom 
element  exhibits  a  fairly  constant  temperature  near  the  centre,  the  top 
element  varies  by  about  25  C  over  the  region  within  4  cm  of  the  centre. 

Thus  the  actual  temperature  differential  within  the  furnace  will  differ  from 
the  applied  value  (Nominally  50  C  during  the  epilayer  growth  runs  reported 
in  reference  1 ,  and  for  the  temperature  probing  measurements  described  in 
this  report)  by  up  to  25°C  over  a  distance  of  just  a  few  cm. 

Since  the  furnace  differential  is  not  constant  even  over  the  central 
region,  then  small  differences  from  run-to-run  in  positioning  the  crucible 
and  substrate  holder  (or  probe)  will  result  in  a  different  temperature 
gradient  within  the  growth  solution.  In  addition,  differences  in  temperature 
gradient  will  occur  over  the  substrate  area,  with  the  side  furthest  from  the 
furnace  being  at  a  higher  temperature  (or,  alternatively,  thought  of  as  under 
the  influence  of  a  reduced  solution  temperature  gradient)  than  the  side  near 
the  centre  of  the  furnace.  It  will  be  recalled  from  Section  3  that  the  grown 
layers  were  in  general  thickest  on  the  side  near  the  centre  of  the  furnace 
(i.e.  near  the  graphite  peg),  and  that  when  etching  occurred  in  Runs  16  to  20 
the  etching  was  greater  on  the  side  away  from  the  graphite  peg.  These  experi 
mental  observations  are  consistant  with  the  observed  reduction  in  temperature 
gradient  and  increase  in  average  furnace  temperature  with  increasing  distance 
from  the  centre  of  the  furnace  elements. 


5.  CONCLUSIONS 
/  AT  \ 

The  vertical  temperature  gradients  (  _ s  )  which  were  measured  within  the 

'  Ax/ 

gallium  melt  do  not  account  for  the  observed  thicknesses  of  the  grown  epitaxial 
layers  if  diffusive  transfer  of  gallium  arsenide  through  the  gallium  melt  is 
considered  as  the  only  process  contributing  to  the  epitaxial  growth.  The  vari¬ 
ability  in  thickness  of  these  layers,  not  only  from  sample  to  sample  but  in  fact 
over  the  area  of  each  individual  sample,  quite  strongly  suggests  that  additional 
processes  are  occurring  within  the  melt. 
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A  plausible  explanation  for  the  observed  situation  is  that  kinetic  processes 
such  as  convection  might  be  taking  place.  The  axial  variations  within  the  melt, 
which  were  found  in  Section  4,  may  set  up  instabilities  which  result  in  cellular 
convection.  Because  there  are  both  horizontal  and  vertical  temperature  variations 
within  the  melt,  and  the  solubility  of  arsenic  in  gallium  is  temperature  dependent, 
the  solute  concentration  will  vary  horizontally  as  well  as  vertically.  The  trans¬ 
fer  of  solute  material  to  the  growth  interface  will  be  increased  by  this  natural 
convection,  thus  increasing  the  growth  rate.  The  fluid  motion  may  create  unstable 
conditions  at  the  interface,  with  any  point  on  the  surface  of  the  growing  layer 
being  subject  to  a  supersaturation  which  varies  with  time.  Such  variations  could 
explain  the  variable  thickness  of  the  epilayer  across  a  given  sample,  and  also 
the  irregular  surface  morphology  which  is  observed. 
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TABLE  1 .  DETERMINATION  OP  THE  SOLUTION  TEMPERATURE  GRADIENT 
USING  THE  SINGLE  THERMOCOUPLE  PROBE 


Day 

Technique 

Solution  (°C/cm) 
temperature  gradient 

1 

lowering  probe 

1.19 

2 

lowering  probe 

0.77 

3 

lowering  probe 

0.53 

raising  probe 

0.74 

4 

lowering  probe 

0.82 

raising  probe 

1.05 

TABLE  2.  SERIES  »A*  MEASUREMENTS  WITH  DIFFERENTIAL  THERMOCOUPLE 
PROBE  (NOMINAL  TEMPERATURE  800°C) 


Furnace 

Absolute  readings 

Relative  readings 

Baseline 

settings  (mV) 

probe  (  yV) 

probe  (  yV) 

drift 

(t  ), 

pi 

(t  ) 

V  p;2 

AT 

P 

Vo 

(  yV/hour) 

V  Ax  / 

4/7/75 

32- 34 

33- 33 

58.9 

48.1 

17.1 

2.05 

+2.5 

5/7/75 

32- 34 

33- 33 

44.3 

8.9 

35.4 

4.25 

-1.3 

•>+  to 
to  to 

1  1 

C\J  to 
to  to 

39.2 

3.8 

35.4 

4.25 

6/7/75 

32- 34 

33- 33 

47.7 

10.4 

37.3 

4.48 

+2.6 

33-33 

32-34 

64.8 

25.0 

39.8 

4.78 

7/7/75 

32- 34 

33- 33 

102.5 

64.3 

38.2 

4.58 

+2.0 

33-33 

32-34 

112.5 

70.0 

42.5 

5.10 

8/7/75 

32-34 

110.1 

36.9 

44.4 

4.29 

5.15 

33-33 

32-34 

117.6 

73.2 

+5.0 

32  mV  setting  corresponds  to  769°C 

33  mV  setting  corresponds  to  793°C 

34  mV  setting  corresponds  to  818°C 


So  that  32-34  represents  a  nominal  furnace  differential  of  50°C 
while  33-33  represents  zero  furnace  temperature  differential. 
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TABLE  3.  SERIES  'B'  MEASUREMENTS  WITH  DIFFERENTIAL 

THERMOCOUPLE  PROBE  (NOMINAL  TEMPERATURE  800°C) 


Furnace 
settings  (mV) 

Absolute  readings 
probe  (  yV) 

Relative  readings 
probe  (  yV) 

Baseline 
drift 
(  yv/houx) 

(Vi 

( V2 

AT 

P 

f AT\  °C/ cm 
^A x  / 

9/7/75 

33-33 

32- 34 

33- 33 

33-33 

32- 34 

33- 33 

361 .5 

367.1 

345 

345 

350 

351.8 

16.5 

16.5 

17.1 

15.3 

1 .99 

1  .99 

2.04 

1  .84 

+1 .8 

TABLE  4.  SERIES  'C'  MEASUREMENTS  WITH  DIFFERENTIAL 

THERMOCOUPLE  PROBE  (NOMINAL  TEMPERATURE  800° C) 


Furnace 
settings  (mV) 

Absolute 

probe 

readings 

(  yv) 

Relative  readings 
probe  (  yV) 

Baseline 
drift 
(  yV/hour) 

(Vi 

(T  )? 

P  2 

AT 

P 

/  AT  \  °C/ cm 
\  Ax  / 

1 5/7/75 

33-33 

32-34 

198.6 

179.7 

18.9 

2.27 

32- 34 

33- 33 

219.0 

199.5 

19.5 

2.34 

+8 

33-33 

32-34 

231 .0 

209.7 

21  .3 

2.56 

1 8/7/75 

Large  drift,  but  average  of  4  changes 
gives 

4.04 

+1 1 

1 9/7/75 

33-33 

32-34 

235.6 

199.0 

36.6 

4.65 

Zero 

33-33 

32-34 

240.3 

195.0 

45.3 

5.43 

32- 34 

33- 33 

227.4 

186.0 

41.4 

4.95 

These  3  runs  on  1 9/7/75  were  done  with  the  probe  in 
different  positions  in  the  melt. 
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TABLE  5.  SERIES  'C'  MEASUREMENTS  WITH  DIFFERENTIAL 

THERMOCOUPLE  PROBE  (NOMINAL  TEMPERATURE  700°C) 


Furnace 
settings  (mV) 

Absolute  readings 
probe  (  y  V ) 

Relative  readings 
probe  (  yV) 

Baseline 
drift 
(  yV/hour) 

(Vi 

(Tp>2 

AT 

P 

/ AT  \  °C/ cm 

^A 7/ 

18/7/75 

29-29 

28-30 

194.0 

159.0 

- - 

35.0 

4.20 

Zero 

TABLE  6.  MEASUREMENTS  IN  HYDROGEN  WITH  DIFFERENTIAL 

THERMOCOUPLE  PROBE  (NOMINAL  TEMPERATURE  800°C) 


Furnace 
settings  (mV) 

Absolute  readings 
probe  (  yV) 

Relative  readings 
probe  (  yV) 

Baseline 
drift 
(  yV/hour) 

(T  ) 
pi 

(Tp>2 

AT 

P 

/AT  \°C/cm 

W  / 

1  9/7/75 

33-33 

32-34 

229.2 

162.6 

66.6 

8.0 

Zero 
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Figure  1 


Figure  1 .  Schematic  diagram  of  the  growth  vessel  used  for  epitaxial 
layer  growth  by  the  Steady  State  method 
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Figure  2 


800°C  for  2  hours 


Top  =  768° 
Bottom  =  817° 


"Average"  value 
K  1 5  p/ hour 


Run  5  Ax  =  4  mm 


/ 

14  10  14  24  24 

27  48 

48 

48  89 

1 

m 

m 

31 

55 

i 

L 

68 

89 


41 


61  34  24  14  20 


Run  10  A  x  =5 


mm 


Run  8  Ax  =2.6  mm 


Run  9  Ax  =  2.5  mm 


|  [68  61 

61 

68 

88 

■1 

15  27 

Run  15  Ax  =  5  mm 


33 

\ 
■'■'30  FAIRLY 

UNIFORM 

28 

29 

~  28  FAIRLY 

UNIFORM 

700°C  for  4  hours 


Top  =  673° 
Bottom  =  721° 


"Average"  value 
~  5  /i  /hour 


Run  11  A  x  =  5  jujq 


Run  14  Ax  =  5  mm 


Figure  2.  Schematic  representation  of  spot  heights  for  epitaxial  layers  grown 
at  "average"  temperatures  of  800°C  and  700°C.  The  left  hand  side 
of  each  diagram  represents  the  edge  of  the  substrate  which  was 
closest  to  the  graphite  peg  used  to  hold  the  source  block  in 
position  in  the  bottom  of  the  growth  vessel 
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Figure  3 


Non-uniform 
growth  (close 
to  graphite 
peg) 


Etching  in  hotter 
region  (distant 
from  graphite  peg) 


Figure  3*  Example  of  non-uniform  epitaxial  growth  and  simultaneous 
etching  of  the  substrate  surface,  caused  by  large 
horizontal  temperature  gradients  in  the  melt  solution 
(Sample  17,  nominally  700°C  for  4  hours) 


Divisions  on  Pen  Recorder 


ERL-0077-TR 
Figures  4  &  5 


Figure  4.  Schematic  of  the  temperature  measuring  system  used  with 
the  single  thermocouple  probe 


Depth  below  surface  of  melt  (mm) 

Figure  5.  Typical  example  of  the  variation  in  results  when  the 
single  thermocouple  probe  is  lowered  into  or  raised 
up  out  of  the  melt  (40  divisions  on  the  recorder 
=  100  nY  from  probe  -  2.4°C) 
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Figures  6  &  7 


pyrolytic  graphite 
probe 

graphite  screw  to 
hold  thermocouple 
in  position 

thermocouple  wires 
encased  in  thin 
quartz  tubing 


Figure  6(a).  Schematic  of  the 

differential  thermo¬ 
couple  probe  assembly 


Figure  6(b).  Enlarged  view  of 

the  differential  probe. 
The  vertical  distance 
between  the  two  chrome 1- 
alumel  junctions  is 
2.0  mm 


Differential  fsd  1  V  fsd 

Thermocouple 


Figure  7 


Technique  used  to  increase  the  measurement  accuracy 
of  the  differential  thermocouple  probe 


Probe  reading  (f/V) 
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Figure  9 


.  Magnitude  of  the  baseline  "zero"  drift  on  one  occasion, 
when  both  top  and  bottom  furnace  temperature  controllers 
were  set  on  33  mV  (793°c)  (average  drift  over  the  linear 
region  shown  is  10.6  juV/Wur) 


Figure  9 
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Figure  1 1 .  Temperature  "overshoots"  produced  in  the  furnace 
by  alteration  of  the  top  or  bottom  furnace 
temperature  controllers 


Thermocouple  Reading  (mV)  (l  mV  approximately  equal  24°C) 


ERL-0077-TR 
Figure  1  2 


Distance  from  centre  of  furnace  elements  (cm) 


Figure  12.  Axial  variation  of  temperature  along  the 

furnace  sections  when  both  controllers  are 
set  to  33  mV  (corresponding  to  793  W 
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